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Abstract 

We apply a boost-invariant similarity renormalization group procedure to a light-front 
Hamiltonian of a scalar field of bare mass fi and interaction term ~ g(j)'^ in 6 dimensions using 
3rd order perturbative expansion in powers of the coupling constant g. The initial Hamilto- 
nian is regulated using momentum dependent factors that approach 1 when a cutoff parameter 
A tends to infinity. The similarity flow of corresponding effective Hamiltonians is integrated 
analytically and two counterterms depending on A are obtained in the initial Hamiltonian: a 
change in /x and a change of g. In addition, the interaction vertex requires a A-independent 
counterterm that contains a boost invariant function of momenta of particles participating in 
the interaction. The resulting effective Hamiltonians contain a running coupling constant that 
exhibits asymptotic freedom. The evolution of the coupling with changing width of effective 
Hamiltonians agrees with results obtained using Feynman diagrams and dimensional regular- 
ization when one identifies the renormalization scale with the width. The effective light-front 
Schrodinger equation is equally valid in a whole class of moving frames of reference including 
the infinite momentum frame. Therefore, the calculation described here provides an interesting 
pattern one can attempt to follow in the case of Hamiltonians applicable in particle physics. 
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1. INTRODUCTION 



Similarity renormalization group procedure for Hamiltonians is a method suggested for 
seeking solutions to QCD and other quantum field theories en bloc including bound states, by 
calculating manageable effective Hamiltonians and solving Schrodinger equations with them in a 
series of successive approximations of increasing accuracy. This article describes an elementary 
application of the similarity method in case of scalar fields, showing details of a calculation of one 
boost-invariant running coupling constant in effective Hamiltonians in third order perturbation 
theory. We briefly review the method, describe the simple example and present conclusions. 

The method originates in the notion of renormalization group as discussed by Wilson [|l| 
and uses the idea of similarity renormalization group procedure for Hamiltonians ^ [^]. 
Similarity enables us to avoid small energy denominators in perturbative evaluation of effec- 
tive Hamiltonians. The evaluation includes finding counterterms and defining renormalized 
dynamics. Small energy denominators could lead to large errors in the counterterms and in 
calculation of effective Hamiltonians. Thus, if not avoided through similarity, the small de- 
nominators would prevent precise theoretical predictions based on the effective Schrodinger 
equations. 

The similarity procedure was invented to sort out complexities of the light-front form of 
Hamiltonian dynamics. This form was distinguished a long time ago by Dirac and more 
recently became a natural candidate for description of hadronic constituents in hard scattering 
processes as well as in spectroscopy using QCD. A recent review article by Brodsky, Pauli 
and Pinsky provides a description of theoretical advances made in light-front formulation of 
various theories mainly before invention or independently of the similarity procedure. Reviews 
by Burkardt, Harindranath and Perry help in understanding the scope of current approaches. 
Recent research in the direction of renormalization of Hamiltonians can be traced through Ref. 
||T0| . The present article is focused on similarity in light-front dynamics. 

Initial studies of quarkonium bound states, which are related to the similarity program 
described in have been performed by Brisudova and Perry, and Brisudova, Perry and Wilson 
TT| following the key observation by Perry [O that second order effective Hamiltonian of QCD 



contains a confining term, which may remain uncanceled in the effective Schrodinger equation 
dynamics. Higher than second order calculations are needed for verifying this hypothesis. 

Since the formal front form of Hamiltonian dynamics is invariant under boosts one hopes 
it can provide a link between the structure of hadrons at rest and in the parton model. That 
such unifying picture is hard to achieve in standard dynamical schemes is best illustrated by 
the fact that despite extensive progress lattice gauge theory does not easily yield desired quark 
and gluon bound state wave functions. |]T3 The light-front approach is still far from achieving 



this goal, too. The present article shows the essence of boost invariance in similarity fl^ but 
the example we describe here is limited to scalar particles. The dynamics of scalars does not 
involve genuine small- a; singularities {x denotes a longitudinal momentum fraction carried by a 
particle in the infinite momentum frame) that appear in gauge theories and the present work 
does not describe known particles. 

The present paper also does not cover the step of solving the effective Schrodinger equation. 
It was shown in an asymptotically free matrix model that one can achieve 10 % accuracy 
in calculating bound state properties using second order effective Hamiltonians but it is not 
known yet if the same accuracy can be reached in any quantum field theory of interest using 
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the boost-invariant similarity approach discussed here. 

Before we proceed to our example, we first disclose the mechanism of preserving boost 
invariance in the similarity approach. The following discussion points out relevant features of 
the procedure using other methods for comparison. 

The similarity renormalization group procedure leads to effective Hamiltonians Hx = Hqx + 
Hjx, whose matrix elements between any two eigenstates of Hqx vanish when the eigenvalues of 
Hqx for these states differ by more than the "width" A. The word "width" is natural because 
the effective Hamiltonian matrix can be viewed as a band of non-zero matrix elements along 
the diagonal and the width of the band depends on A. Another reason for the word "width" 
is that the band structure is ensured by similarity form factors in interaction vertices. The 
form factors are functions of energy transfers, they are peaked around zero and their half-width 
depends on A. The notion of the Hamiltonian width is key to our method of preserving boost 
invariance. 



We shall take advantage of Wegner's equation ||T6[ to indicate where the boost invariance 
can enter the similarity procedure. Wegner invented a flow equation for diagonalization of 
Hamiltonian matrices in solid state physics ||16[ that is beautifully simple and can be 



adapted to the similarity renormalization scheme [|T8| [|T^. Wegner's equation for Hamiltonian 
matrices can be written as 

^^Hx = - [[Hx,Hox],Hx] , (1.1) 

where s = A~^. Initial condition should be provided at s = 0, corresponding to A = oo, and the 
initial Hamiltonian is denoted by Hoo. Wegner's commutator [Hx, Hqx] generates the similarity 
transformation. One discovers a gaussian similarity factor by solving Eq. (1.1) for the Hjx 
matrix elements between eigenstates of matrix Hqx keeping on the right-hand side only those 
terms that are linear in Hjx and neglecting higher order terms. 



Wegner's equation preserves necessary cluster decomposition properties [|l^]. The commu- 
tator structure of Wegner's equation implies that no disconnected interactions are generated 
by the transformation. Wegner's transformation depends only on differences of energies and 
spectator energies always drop out from the differences. However, Wegner's generator is not 
boost invariant. 

Seeking a boost invariant approach, one can apply Eq. (1.1) to light-front Hamiltonian 
matrices. The matrices are of the form H = (P-*-^ + M^)/P~^, where the mass matrix contains 
interactions, M| = Mqx + Mjx, and P denotes total momenta of states with which one evaluates 
the matrix elements. Since the masses do not change the momenta, one can rewrite Eq. (1.1) 
in terms of the mass matrix elements as follows 

P^'j-^Ml = - [[Ml M'qx], Ml] . (1.2) 

We see that a rescaling of the flow parameter s with a momentum eigenvalue gives a flow 
equation for the mass squared matrix elements. The latter should be independent of the eigen- 
values P+ and P^. This feature of masses seems to suggest a boost invariant renormalization 
procedure based on replacing P^"^ / s by a new P+-independent flow parameter A^. But such 
substitution breaks connection with Wegner's equation for Hamiltonian matrices and leads 
to violation of cluster decomposition properties, since the mass is not an additive quantity. 
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Namely, the mass depends on the relative motion of particles and effective interactions become 
dependent on spectators momentum. 



Despite this drawback, Allen and Perry [^] succeeded in demonstration that one can define 
and calculate a running coupling constant using Eq. (1.2) for matrix elements replacing 
P^'^/s with in massless 0^ theory in 6 dimensions. The authors demand that the flow from 
some A to a fraction of A reproduces the same matrix elements structure in which only some 
parameters change. This condition is implemented using transverse locality and it allows for 
bypassing the step of finding initial conditions (counterterms) needed for differential equations, 
by introducing a running coupling constant. A question arises because studies of asymptotically 
free matrix models |jl5| show that effective Hamiltonian matrices that are suitable for 



bound state calculations may significantly deviate from a self- replicating (fixed point) flow 
with one coupling constant. But in the case of many couplings helpful conditions may be 
provided by coupling coherence which may work in the approach of Eq. (1.2) with a P~^- 



independent width parameter, assuming that the cluster property breaking does not complicate 
the coherence. Kylin, Allen and Perry |^ extended the approach of Allen and Perry to massive 
particles. 

The similarity renormalization group procedure for Hamiltonians is more flexible in defining 
its generator than Eqs. (1.1) and (1.2) might suggest. This flexibility is used to preserve boost 
symmetry and cluster properties simultaneously. Note that the boost invariance we talk about 
means invariance under boosts along the front and boosts corresponding to rotations about 
transverse axes in the infinite momentum frame. These boosts are sufficient for reaching states 
in all possible uniform motions. 



The desired boost-invariant operator formulation of similarity is based on a transfor- 
mation U\ that changes creation and annihilation operators of bare particles (denoted by goo, 
since they correspond to A = oo) into creation and annihilation operators of effective particles 
corresponding to the width A (denoted by ^a)- Namely, 

qx = Ux Qoo Ul . (1.3a) 

lAx is secured to be unitary by construction. Hamiltonian operators of all widths are assumed 
equal and when they are expressed in terms of different creation and annihilation operators, 
the coefficient functions change. We have Hx{qx) = -f^oo(g'oo)- Assuming that Hamiltonians 
calculable in perturbation theory contain only finite products of creation and annihilation 
operators and applying the transformation Ux, one obtains Hx = Hx{qoo) = Wa-^oo(?oo)Wa. 
This relation means that the operator Tix has the same coefficient functions in front of products 
of qoo as the effective Hamiltonian Hx has in front of the unitarily equivalent products of gA- 
Differentiating Tix one obtains 



d 
J) 

where the generator Tx is related to Ux by 



— nx = - [Ta, Ha] , (1.36) 



Tx = ui^Ux . (1.3c) 

The script letters are introduced in Eqs. (1.3a-c) to indicate that the operators can be con- 
veniently thought about as expanded into sums of products of operators goo- The latter are 
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independent of A and are not differentiated in Eqs. (1.3b-c). In other words, Eqs. (1.3b-c) de- 
scribe only the flow of coefficients in front of the creation and annihilation operators. Effective 
Hamiltonians are obtained from Hx using Hx{qx) = Ux^x^x- 

Note that the operator Ux is analogous to the Melosh operator . The analogy stems from 
that the effective particles corresponding to a small Hamiltonian width A can be associated with 
constituent quarks or gluons from hadronic rest frame spectroscopy while the bare particles with 
A ^ oo can be associated with current quarks and gluons from the infinite momentum frame. 
The boost-invariant similarity renormalization group procedure for light-front Hamiltonians 
|l^ makes the notions of current and constituent quarks and gluons boost invariant. The 



procedure provides a calculable dynamical content to the transformation that connects the 
current and constituent particles independently of the reference frame one works in. This may 
sound confusing if Melosh transformation is associated with changing a frame of reference. But 
the light-front dynamics is invariant with respect to boosts in perturbative calculations and this 
property allows for translating a Melosh-like transformation to one frame of reference, where it 
simply relates bare particles and complex effective particles in one theory in the same frame. 

Boost invariance is guaranteed through the definition of the generator Tx- Details will be 
reviewed in the next Section. The effective Hamiltonian Hx is given by a diagonal proximum of 
certain operator Gx, i-e. Hx = Fx[Gx], while the generator Tx is related to a diagonal remotum 
of the same operator Gx- Following |]l4l, we introduce Qx = uIGxUx and use 



[Tx^mx] = ^{l-Fx)[gx] . (1.4) 

The commutator structure guarantees that effective Hamiltonians contain only connected in- 
teractions. 

Now, the diagonal proximum can be defined in a boost-invariant way. The point is that 
the operation Fx multiplies every interaction term in an effective Hamiltonian by a form factor 
which is a function of boost invariant combinations of momenta labeling creation and annihi- 
lation operators in that interaction term. Such combinations can include differences or sums 
of invariant masses. The resulting generator of the similarity transformation does not depend 



on the total momentum of any state, contrary to Eq. (1.2). According to [|T4|, a necessary 
rescaling of the flow parameter, analogous to rescaling provided by the total P+ in Eq. (1.2), 
is provided for each Hamiltonian term separately using a parent three-momentum for that 
term. The parent momentum is defined as half of the sum of all momenta labeling creation 
and annihilation operators in the term in question. All light-front Hamiltonians we consider 
are sums of terms which contain integrals over the parent momenta. Usually, when a term acts 
on some state, a parent momentum equals only some fraction of the total momentum of that 
state. This fraction is defined by the single interaction act and lies in a range of values allowed 
by momentum conservation in the effective Schrodinger dynamics. Thus, effective interactions 
are independent of the total momenta of states they act on and no dependence on spectators 
is generated. Therefore, no cluster property is violated in defining the similarity generator 
through Eq. (1.4). A considerable freedom is still left in choosing details of Fx- 

In summary, the flexibility available in defining the similarity transformation generator can 
be used to obtain a boost-invariant band-diagonal structure of effective Hamiltonians preserving 
cluster decomposition properties, through a suitable choice of the similarity form factors Fx- 
This aspect of the procedure is the subject of the present article, quite independently of which 
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theory or singularity is addressed. Naturally, the effective theory is easiest to make boost 
invariant if the initial Hamiltonian is regulated in a boost invariant way. For in that case the 
counterterms do not have to correct the boost invariance breaking which a frame dependent 
regularization would introduce. A suitable class of regularizations will be described in the next 
Section. 

This paper is organized as follows. Section 2 reviews necessary elements of Ref. [|T3| set 



up for a third order calculation of the running coupling constant in Hamiltonians of scalar 
(f)^ theory in 6 dimensions. The calculation is based on a plain power series expansion in the 
coupling constant. Counterterms are derived in Section 3 and the effective coupling constant 
flow is calculated in Section 4. Section 5 concludes the paper. Appendix presents formulae for 
third order effective vertex function with arbitrary momenta and masses. 

2. EFFECTIVE HAMILTONIANS 

We focus our attention on a perturbative derivation of effective Hamiltonians in scalar 
quantum field theory with an interaction term ~ 0^ in 6 dimensions. The theory is known to 
be asymptotically free. The coupling constant in the initial regularized Hamiltonian is assumed 
to be infinitesimally small. Our procedure for calculating effective Hamiltonians closely follows 
Ref. 0. 



2.1 Initial Hamiltonian 

The classical Lagrangian of the scalar field theory in question is 



£=^(9^09>-/xV)-|0' • (2.1) 



The corresponding light-front Hamiltonian reads |0 ||25| (x^ = ± a; 



-'-J./ n±2 I ,,2\i I 9 ,3 



H = I dx d'^x^ 

x+=0 

A quantum field (f){x) for = is expanded in its Fourier components, 



(2.2) 



0(x)U+=o= / [A;](afce-*^- + 4e*'=-)U+=o , (2.3) 
where the abbreviated notation for the integral means 



The small parameter 5^ limits the longitudinal momenta. A;"'", from below. The creation and 
annihilation operators satisfy commutation relations 

K, aj] = 2fc+(27r)55(fc+ - q+)5\k^ - q^) , (2.5) 

together with [a^, Og] = [a\^ aj] = 0. Substituting Eq. (2.3) in (2.2) one obtains the following 
unique structure of the light-front Hamiltonian 
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H = jjk] ^ ^^^^ + lj^^kP<iW^f^'ik + p-q){aiala, + a\a,a,) . (2.6) 

It is unique in comparison to Hamiltonians in other forms of dynamics because there are no 
terms that contain only creation or only annihilation operators in Eq. (2.6). This feature 
is related to the problem of ground state formation, since H\0) = 0, where ak\0) = 0, and 
only alalal and terms in H could alter this feature, but they are absent due to the cutoff 
S~^. Readers interested in non-perturbative aspects of Hamiltonian dynamics in the front form 



should consult p6[. However, no problems with modes of A; = arise in perturbation theory 
for Hamiltonians in the present paper. 

The initial Hamiltonian for the similarity renormalization group procedure is obtained from 
Eq. (2.6) by introducing an additional ultraviolet regularization factor (to be described below) 
and taking the limit 5+^0. The initial Hamiltonian is denoted by i^oo? since it corresponds 
to the initial width A = cxd in the similarity flow, as discussed in the Introduction. Creation 
and annihilation operators that appear in H^o were denoted in the Introduction by goo- 

i^oo contains a regularization factor denoted by r^. A stands for an ultraviolet cutoff 
parameter, ta — * 1 when A ^ oo. In asymptotically free theories, A can be sent to infinity 
when the renormalization process passes the stage of calculating counterterms and deriving 
effective Hamiltonians of finite width A. Still, the initial Hamiltonian contains counterterms 
which depend on ta and are denoted by Xa- Besides diverging A-dependent terms, Xa contains 
also finite parts that remove A-independent regularization effects caused by the factor ta. Thus, 
we have 

fmf^^^ + f^^ t 
Hoc = [k\ 7- a^ofc^oofe + 



+ I y [kihh] 2{2-K)^6^{ki + ^2 - ks) (aLfcioLfca^oofcs + aLfcsOoofcaOoofci) ^a + Xa ■ 

(2.7) 

The parameter is set equal and this is why it is not indicated, in distinction from Eq. (2.6) 
where it was kept larger than 0. 

The regularization factor ta has a simple form which results from the following steps (ap- 
plicable to the interaction term written out explicitly in Eq. (2.7) as well as to all counterterms 
Xa derivable in perturbation theory). For a term containing a product of u creation and v 
annihilation operators we define a parent momentum, denoted by Pj^, which equals half of 
the sum of momenta labeling all the operators. For each momentum label fcj, with i running 
through u + v values, we introduce Xi = kf / P^^ and = kf- — XiP^^. The regularization factor 
is defined as 

_r). 

rA = n^xp^^ , (2.8) 

1=1 ^ 

where rji = r]{xi) and ?7 is a useful function of its argument. One natural choice is r]{x) = 1, 
for it is simple. Another choice is ri{x) = l/x, a natural one because it appears in invariant 
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masses. We shall assume the function rj lies between the two choices but otherwise it will be 
left unspecified. Leaving t] unspecified will help us identify finite regularization effects. 
Consequently, the regularization factor in the second term in Eq. (2.7) is 



A2 



rA = exp Vr'"^'' ' (2-9) 

where Xi = kf /k^ and — ~ Xik^. 

The initial liamiltonian of the similarity renormalization group procedure in 0^ theory in 6 
dimensions is given by Eq. (2.7). In order to derive the third order running coupling constant 
in effective Hamiltonians, we have to calculate up to the terms order g^. The counterterms 
are calculated order by order along the evaluation of renormalization group flow for effective 
Hamiltonians. 



2.2 Similarity Flow of Hamiltonians 

The effective Hamiltonians are written as 

Hx = Fx[Gx] , (2.10) 

where 

Gx = UxQxUl . (2.11) 

The operator Qx is divided into two parts, Qx = Go + Qix, where Qq = Gooio = 0). Qix satisfies 
the following differential equation as a consequence of Eqs. (1.3a)- (1.4) (see Ref. pi| ) 



(2.12) 



We dropped the subscript A on the right-hand side for clarity. / denotes the similarity form 
factor introduced by Fx and the curly bracket with the subscript Qo denotes a solution for Tx 
resulting from Eq. (1.4). We will omit the subscript Qo from now on and the Reader should 
remember that a curly bracket implies an energy denominator for every term in the bracket, 
i.e. a factor equal to inverse of the eigenvalue of Qq corresponding to momentum labels of all 
annihilation operators in a term minus eigenvalue of Qq corresponding to all creation operators 
in the term. 

The similarity factor is defined for any operator of the form described above Eq. (2.8) as 

fxiu,v) = exp [-{Ml- Mir /X""] . (2.13) 

The script notation for invariant masses means M^ = {ki + ... + fc„)^, where the minus compo- 
nents of the momentum four-vectors are given by k~ = {k^"^ + ^'^)/kf for « = 1, u. 

Let us denote differentiation with respect to A by a prime and expand the effective interaction 
in powers of the coupling constant g as 

oo 

<3i = E ^" ' (2-14) 

n=l 
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where r„ denotes the sum of all terms order g" in Qx- Equation (2.12) implies 

n-l 

< = T.[n,{a-f)rn-ky] , (2.15) 
fc=l 

and the first three terms in the expansion satisfy the following equations 

t[ = , (2.16a) 

< = [{f'njjn] , (2.166) 

^3 = [fn, {-/V2 + (1 - /)r^}] + [/r2, {-f'n}] . (2.16c) 

Equation (2.16a) implies that ti is independent of A and equals the second term in the 
initial Hamiltonian from Eq. (2.7). In other words, r^i = Tooi- The corresponding effec- 
tive Hamiltonian interaction term is obtained by multiplying the integrand in Eq. (2.7) by 
/a(12,3) = exp { — [(fci + ^2)^ — /^^]^/A^} and transforming operators Ooofc into axk using Ux 
after the evaluation process for Qx is carried out to the desired order. 

This last step is unusual in the sense that the operator Ux depends on the regularization 
(in scalar theory, the dependence is reduced to Ux being a functional of ta) but, at the same 
time it is unitary and thus, it transforms finite products of creation and annihilation operators 
by effectively replacing everywhere Ooo by ax and no other trace of Ux is left in the effective 
Hamiltonian. If Qx is found to order g^, leaving terms o{g"'~^^) still undetermined, then Tx 
is determined up to terms order g"' from Eq. (1.4). 7^ is antihermitian order by order and 
determines a unitary Ux to order g"-, denoted U^\ Therefore, 

^1"^ n«Ln«oo,wf^' = Il4^I[ax, + oig-^') , (2.17) 

and the regularization dependence of Ux does not show up in Hx to order g"- once Qx is renor- 
malized to order g"'. In successive orders, by construction, the counterterms in Qoo preserve 
unitarity of Ux- A perturbative proof of renormalizability for effective Hamiltonians would 
require demonstration that there exists a set of counterterms that remove regularization de- 
pendence from finite momentum parts of Qx to all orders of perturbation theory when A (yo 
[c.f. 0) and that the resulting theory predicts covariant results. The present article does not 
demonstrate such set exists in the case of 0^ theory. Also, our calculation is limited to terms 
order g, g"^ and g^. To verify that an effective Hamiltonian containing a running coupling con- 
stant order g^ may produce covariant results for scattering processes, the present calculation 
must be extended to 4th order. 

The transformation connecting two effective Hamiltonians with different finite widths Ai 
and A2 is given by Ux^Ul^. One can easily see that the latter is free from dependence on ta 
once Qx is made independent of r^. It is sufficient to observe that for an infinitesimal difference 
between Ai and A2 the effective transformation Ux-JJx2 given by the similarity generator 
that is expressed in terms of the rA-independent Qx- Integrating the infinitesimal changes one 
obtains the same conclusion for finite changes of the width. 
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Evaluation of T2 and T3 involves calculation of two counterterms. Both T2 and are more 
complicated in structure than ti and their evaluation requires new notation. Namely, 



n = + ai2 , (2.18a) 

r2 = + /931 + /?13 + /522 , (2.186) 

T"3 = 721 +712 + 741 +714 + 732 +723 ■ (2.18c) 

Each term contains products of creation and annihilation operators with fixed numbers of the 
operators in a product. The first subscript indicates the number of creation operators, a}^, 
and the second subscript denotes the number of annihilation operators, aoo- For all the terms. 



Equations (2.16a-b) imply for second order terms the following relations [|T4 



P31 = 2/2[a2ia2i]3i , (2.19a) 
= 2/2[ai2ai2]i3 , (2.196) 
f3[^ = 2/2[ai2a2i]ii , (2.19c) 

P22 = /2[a2iai2 + 4ai2a2i]22 . {2.19d) 

The brackets mean replacement of products a^a] by commutators [a^, aj]. The remaining prod- 
ucts of operators contain as many creation and annihilation operators as indicated by the 
bracket subscripts, in the normal order, according to the same convention as in Eqs. (2.18a- 
c). The factor /2 depends on the momenta labeling creation and annihilation operators. In 
Eqs. (2.19a-d), the brackets with operators involve integrals over three-dimensional momentum 
labels of all creation and annihilation operators and over loop momenta in loops that result 
from contractions (in second order here only Eq. (2.19c) for f3[^ contains a loop integral). The 
factor /2 is understood to appear under the integrals. Symbolically, its structure appears in 
Eq. (2.16b) and reads /2 = {/'}/ — /{/'}• The negative sign results from the commutator 
in Eq. (2.16b) that guarantees that only connected terms appear in the effective interactions. 
This is a general property of the similarity procedure. 

The factor /2 is the only factor depending on A on the right-hand side of Eqs. (2.19a-d). 
Therefore, solutions are 

/3a3i = 2J'2A [a2i«2i]3i , (2.20a) 
Pm3 = 2J-2A[ai2ai2]i3 , (2.206) 

Pxil = 2J^2A [ai2tt2l]ll + Pool! , (2.20c) 
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P\2i = ^2A [oi2iOii2 + 4q;i2Q;2i]22 , (2.20d) 
where Poou is a counterterm whose structure is shown in the next Section to be 

/3ooii = / [k] ^-^al,a^, . (2.21) 

The only new element of solutions (2.20a-d) which requires explanation is the factor J^2\ — 
f2- It is given by the following expression, which is a consequence of Eq. (2.13), 

^2A(a,6,c) = 'la?Xhc^ [hMh{h,c)-l\ . (2.22) 

Arguments a, h and c denote three successive momentum configurations appearing in brackets in 
Eqs. (2.19a-d) in the order from the left to right, i.e. a denotes the momenta labeling creation 
operators in the brackets, c denotes the momenta labeling annihilation operators and b denotes 
the intermediate configuration, which includes the momenta labeling operators contracted in 
the brackets and momenta labeling creation operators originating from the interaction con- 
necting configuration h with c and momenta labeling annihilation operators originating in the 
interaction connecting configuration a with h. 

P^^ denotes parent momentum for the whole connected interaction sequence between mo- 
mentum configurations u and v (in the second order case here the sequence reduces to a single 
interaction vertex order g^, but the definition of remains valid in higher order cases later). 
The symbols ha and he denote differences of invariant masses, as explained below. We use 
abbreviations ha? = {haY etc. 

For any two momentum configurations u and v 

uv^Ml,-Ml^ , (2.23a) 

where 



•' '■uv 



E k^ (2-236) 

and u{v) denotes those momenta from the configuration u that are involved in interactions 
acting between the configurations u and v. Minus components of all momenta are given by 
A;" = (k^"^ + n'^)/k^ . As an example of this notation, the similarity form factor from Eq. (2.13) 
reads 

2 

— VZt 

fx{u,v) = exp^^ . (2.23c) 

Equation (2.16c) for third order interactions implies the following result for terms relevant 
to the running coupling constant evaluation, 

721 = /3[8[ Q;i2CK2iQJ2i]2i ~I~ 4 [[q;i2Q;2i]q;2i]2i + 2 [q;2i[q;i2Q;2i]]2i ] + 

- 2{/'} [/?ooll«2l]21 + {/'} [a2l/?ooll]21 ■ 
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The factor /a has the structure /a = [/{(I - f)J^2}' + {/'}/-^2] 
Integration of both sides of Eq. (2.24) gives 



(2.24) 

- [{(l-/)J-2}'/ + /J-2{/'}]. 



7A21 = ^3A [8 [ai2a2ia2i]2i + 4 [[Q;i2a2i]a2i]2i + 2 [a2i[tti2Q;2i]]2i] + 



+ 2{1 - fx} [/3oolia2l]21 - {1 - fx} [a2l/5ooll]21 + 7oo21 , 

(2.25) 

where 70021 denotes the third order counterterm, to be calculated in the next Section. A new 
element of Eq. (2.25) is the factor J-'^x = J^f?,- It appears in front of operator brackets that 
involve four successive momentum configurations denoted from the left to right by a, b, c and 
d. The brackets contain one loop integral. The configuration a has two momenta, ki and /c2, 
while the configuration d only one, ks [c.f. Eq. (2.7)]. We have 



J^3X {a,b,c,d) = F3 (a, 6, c, ci) + F3 (d, c, 6, a) 



(2.26a) 



where 



^, , P+cb + P+cd 

F3{a,b,c,d)^ 



dp' + cd? 



{P^,bd + P+6a) 



fabfbcfcdfbd 1 



fabfl 



bd 



+ Pb^d- 



bc^ + of 



db 



fabfbcfcd 1 

a62 + + ~ + bd^ + cd^ + bd"^ 



ab^ + bc^ + cd^ + bd^ ab^ + bd^ 

fabfbcfcdfbd 1 



+ 



(2.266) 



and fab is an abbreviated notation for fx{a,b) = exp [— a6^/A'']; c.f. Eq. (2.23c) above. 
3. COUNTERTERMS THROUGH ORDER 

The two counterterms, /Sooii from Eq. (2.20c) and 7oo2i from Eq. (2.25), are calculated 
using the equations they appear in. The counterterms are determined by the condition that 
those equations are independent of ta when A — > oo for arbitrary finite values of A and particle 
momenta. We first describe calculation of Poou and then 70021 ■ 



3.1 Mass counterterm 

Equation (2.20c) implies 



All 



where 



[k] 



k+ 



^oofc^oofe 



(3.1) 
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*^5-^''^+(f)^^I^/''*''^^^5^[AW..Vl]r., . (3,2a) 

The script Ai denotes invariant mass, Ai^ = (k-*-^ + — x), and the regularization factor 

comes out to be 

TAfs = exp[-2[r]{x) +r]{l-x)]K^yA^} . (3.26) 

in Eq. (3.2a) is the counterterm contribution. The counterterm is of the form given in 
Eq. (2.21) since the integral in the second term in Eq. (3.2a) is independent of the momentum 
k and depends on regularization. The result of integration depends on the cutoff parameter 
A and function rj. A finite, regularization dependent part of the result remains undetermined. 
Therefore, we have to adjust its value to data. 

Without loss of generality, we assume that some gedanken data requires the mass squared 
parameter in effective Hamiltonian with A = Aq to be equal /x^ + (J/ig- Hence, 



2 J 2(27r)5 Jo x{l-x 



dx 



Tap 



(3.3) 



Integration gives two diverging terms, one proportional to and another one proportional 
to log A. The remaining finite part depends on our choice of the function rj. Evaluating the 
integral for ri{x) = 1/x one obtains 



2 

oo 



(4vr) 



1 

24^ 



/i ^log — 

D /i 



(3.4) 



where has a finite limit when A ^ oo. The logarithmically divergent part is independent of 
the function rj and agrees with results for the Lagrangian mass squared counterterm obtained 
using Feynman diagrams and dimensional regularization ||2^ ||28| in the following sense: when 
one changes A to A' the logarithmic part of the counterterm changes with A as the mass 



squared changes as a function of the renormalization scale in Eq. (7.1.22) in |28 



The resulting mass squared term in the effective Hamiltonian can be written in the limit 
A — > cx) as 



2 2 I c- 2 



2 , e 2 , f9\^ 1 dx /• 4 , 2 



2 J 2(27r)5 7o x{l - x) J M"^ - ii^ 



d\' 



(3.5a) 

The above result is particularly simple for /i = and in that case it reads ((5;Uq is proportional 
to g-") 
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Logarithmic dependence on A arises for fi > 0. The value of 6fiQ could be found, for example, 
by solving a single physical meson eigenvalue problem, expressing the physical meson mass 
in terms of S/i^ and adjusting the latter to obtain the gedanken experimental mass value for 
mesons. 

3.2 Coupling constant counterterm 

The couphng constant counterterm 70021 is evaluated from Eq. (2.25). The mass squared 
counterterms cancel quadratic and part of logarithmic divergences so that the vertex countert- 
erm is needed to remove only logarithmically divergent and finite regularization effects. 

The interaction term 7^21 has the form 

7A21 = j [kik2k3\ 2{27rfS^{ki + k2 - ka) [-fx{ki, /c2, k^) + 700(^^1, k2, ks)] alofci^lofea^oo^^s '^a ■ 

(3.6) 

The counterterm function 700(^1, k2, k^) originates from 70021 in Eq. (2.25). Since the similarity 
renormalization group procedure preserves canonical symmetries of light-front Hamiltonians, 
the functions 7a(^i, k2, k^) and '~ioo{ki, /c2, k^) depend only on variables X\ and k^2-i introduced 
in Eq. (2.9). 

The entire regularization dependence of 7a(^i,^2,^3) = 7a(3^1)'^i2) in the limit A — > 00 is 
contained in 

7A(a;i,fi;^2)|rA = 7A(a;i,0^)|^^ = (^-j x 

/■I dx f°° 2 , 2 ^ ~ ^1 ( —c„K^\ 1 
/ r / n an 8 — -exp — \ + [xi X2) + 

JxiX{l-x){x-Xi)Jo XX2M'^ ^\ J 

(3.7a) 

where 

= r]{x) +r]{l-x) + {r]{xi/x) + -q^x - Xi)/x]} {xi/x)^ + r][{x-Xi)/x2] + r][{l - x) / X2] (3.76) 

and 

dr, = 2 [77(0;) + 77(1-0;)] . (3.7c) 
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The first term in Eq. (3.7a), symmetrized in xi and X2, originates from the first 8 terms in 
the long bracket in Eq. (2.25). The second term originates from the next 4 + 2 terms in 
the long bracket and 2 + 1 mass counterterm terms in Eq. (2.25). The regularization effects 
are independent of k^j- The apparent singularity at — > for = is irrelevant to the 
regularization dependence and appears here only because we do not display similarity factors 
that remove the singularity. The full expression for 'y\{xi, k^j) is given in Appendix. 
Dropping all parts that are independent of regularization, Eq. (3.7a) gives 



-fxiXi,0^)\r^ 



A 

3 log 4 

/i 



,27 (47r)3 
"i 1-x 



X 



/ dx logc„ + (xi X2) + 3 / dx x(l — x)\ogdn 

Jxi X2 J Jo 



(3. 



The bare boson mass, fi, can be replaced by a finite constant of the same dimension in the case 
of massless bosons. 

Equation (3.8) says that the diverging regularization dependence of the interaction vertex, 
i.e. the term proportional to log A, is independent of the particle momenta and one can 
remove the divergence by merely changing the initial coupling constant g in Eq. (2.7). Thus, 
no diverging x-dependent counterterms are required - different situation than in However, 
it is visible that the vertex contains a finite regularization dependent part that is a function of 
Xi. The function depends on our choice for 77. One could completely subtract the ?7-dependence 
for arbitrary choices of rj by defining a counterterm that contains a negative of the 77-dependent 
part of the right-hand side of Eq. (3.8). For example, if 77 = 1 one has = 4 + 2{xi/x)'^ and 
drf = 4. Integration produces a smooth and slowly varying function of xi, which needs to be 
subtracted. 

Since the whole regularization effect is independent of A and k^2! can be completely 
removed from 7a(xi, K;j'"2) by subtracting ■yxo{xi,0'^), where Aq is chosen arbitrarily. However, 
one has to add back the finite regularization independent part of 7Ao(a;i, O"*-), which will be 
denoted below by 7o(a;i). The function 7o(a;i) is necessary to recover Poincare symmetry of 
observables. Regularization spoils Poincare symmetry. The symmetry may be restored once 
counterterms remove regularization effects, but one is not allowed to change terms independent 
of regularization. Therefore, the function 70 (xi) must be reinserted. This function is not altered 
when A changes and could be considered marginal in analogy with usual renormalization group 
analysis. 

Although one can isolate finite 77-dependent functions of Xi in the effective vertex, for 
the particular choice of the regularization factors we adopt, the ultimate adjustment of the 
function 70 (xi) has to be delayed until 4th order calculations are completed. For there exists 
in (j)^ theory no 3rd order scattering amplitude one could use to find out what function 7o(a;i) 
renders Poincare symmetry of scattering observables with our choice of ta in Eq. (2.7). It will 
be interesting to see if the explicit dependence on rj isolates the same function as required by 
the symmetry. However, it should be pointed out that the function does not influence the way 
the 3rd order running coupling constant in effective Hamiltonians depends on A. 
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The counterterm function 700(^1, ^2, ^3) = 7oo(a;i), which removes the regularization depen- 
dence from the effective vertex reads 



7oo(a;i) = -7Ao(a;i,0^) + 7o(a;i) ■ (3.9a) 

It can be used to define a new regularization dependent couphng constant qa in the initial 
Hamiltonian in Eq. (2.7). We select a convenient value of xi = xq and obtain 



9A 

2 



I + 7oo(a;o) 



I - 7Ao(a;o,0^) +70 



(3.96) 



where 70 = 70(3^0)- Using Eq. (3.8), we see that the initial coupling g is replaced by the new 
A-dependent quantity 



9a = 9 



9 



4(47r)' 



log 



A 

mo 



(3.9c) 



with certain constant mo- Thus, the theory exhibits asymptotic freedom in 3rd order terms. 
Our result agrees with literature, say Eq. (7.1.26) from [Q, in the sense that when we change 
A, the change required in the coupling constant in the initial Hamiltonian for obtaining A- 
independent effective Hamiltonians matches the change implied by Feynman diagrams and 
dimensional regularization. Comparison with Feynman calculus will be farther discussed below. 

Having derived the structure of counterterms we can proceed to evaluation of the finite 
similarity flow of effective Hamiltonians towards small widths A. 



4. RUNNING COUPLING CONSTANT 

Our procedure for evaluating the running coupling constant in Hamiltonians follows theory 
from and using . The procedure has been outlined in a matrix model example in [|T5| 



Here, we use particle creation and annihilation operators in a boost invariant way, instead of 
matrix elements. 

The initial Hamiltonian interaction vertex depends on regularization through the factor 
ta and corresponding counterterm. The bare coupling constant g is replaced by the coupling 
constant g/\ according to Eqs. (3.9a-c), i.e. 



9a = 9 -"2 7Ao(a;o,0-^) - 7o 



(4.1) 



The factor 2 is needed because g/2 appears in the Hamiltonian. Both 7Ao(a^o) and 70 are 
proportional to g^. Inverting the series (4.1) one can express g in terms of g^. 

Evaluation of 7A21 in Eq. (3.6) leads now to a finite expression, which has a limit when 
A —>■ 00. In that limit, the effective Hamiltonian interaction term takes the form 



H\2l 

where 



[hk^h] 2{2'Kf5\k^ + k2- ks) fx[{k, + k2)\ J^i] Vx{x^ yi^u) (^Xki(^xk2^\ks , (4.2a 



9A 



Vx{xi,K^2) = ^ + lx{xui^i2) +lxoi^o,0 )-7Ao(a;i,0 )] + 7o(a;i) - 7o 



(4.26) 
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is the effective vertex function and fx is the similarity vertex form factor, ^yxlxi, K12) and the 
remaining terms in Vx{xi, K12) are proportional to g^. Our calculation is done to this order of 
accuracy only, so that = g\ and Eq. (4.2b) is considered to be an expansion in powers of 
Qa- 

We define the running coupling constant as the value of 2Vx{xi, ^^2) at a chosen config- 
uration of momentum variables, denoted by {xio, k,^2o) specified later. In other words, 
gx = 2VA(a;io, /^i2o)- This is a natural definition, analogous to the standard Thomson limit in 
the case of electron charge in QED. We have 

9\ = 9a + 2[7A(a:io, ^.^20) + Ixoi^o, 0"^) - 7Ao(a;io, 0"^) + 7o(a;io) - 7o] ■ (4.3a) 
It is natural to use Xq = Xiq. Then, 

gx = Qa + 27a(xio,k]^2o) ■ (4-3&) 

This equation demonstrates that the effective coupling constant gx depends on the value of the 
finite function 70(2^1) in the counterterm at one point, the same as the one used to define g^- 
Suppose that for the chosen value of A = Aq, the running coupling constant should have the 
value gxo — go, determined from comparison with experiment. Then, 

9o = 9a + 27Ao(a;io,K^2o) + o(^a) > (4-4) 
where g in terms order g^ is replaced by g/^. Inverting this series expansion we obtain 

gA = ^0 - 27ao(xio,k^2o) + 0(^0) > (4-5) 

where in terms order g\ we have g^ replaced by go. Relation (4.5) can be inserted into Eq. 
(4.3b) to yield 

gx = go + 2 [7a(xio,«;^2o) -7Ao(2;io,«;m)] + ^(fifo) ■ (4-6) 

This relation is free from dependence on the finite function 7o(a;i) in the counterterm. 
The vertex function in the effective interaction in Eqs. (4.2a-b) is equal 



Vx{xi,Ki2) = Y+7A(a:;i,«^2)-7Ao(2^io,«^2o)+7Ao(2^io,0^)-7Ao(a;i,0^) + 7o(xi) -70 + 0(^0) > 

(4.7a) 

where in terms order g^ the initial g is replaced by go. The difference between g and go is of 
order g^ and it is included in terms o(gfo)- Equation (4.7a) gives us the effective vertex function 
for width A expanded in powers of the effective coupling constant go corresponding to width 
Aq. Written as a power series in gx, the effective vertex function reads 



Vx{xi,Ki2) = y + 7A(a;i,«^2)-7A(a;io,«m) + 7Ao(a^io,0^)-7Ao(a;i,0^)+7o(a;i)-7o + o(^f) , 

(4.76) 

where in terms order g^ one replaces g by gx. Clearly, Eq. (4.7b) reproduces the relation 
gx = 2Vx{xio,K,i2o)- 
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It remains to calculate the dependence of gx on X, which requires a choice of the momentum 
configuration (xio, t^uo)- The particular choice we will adopt is suitable for massless bare bosons, 
i.e. for /i = 0. Our choice would require a change for /i > 0, to avoid vanishing of similarity 
factors when /i^/xio tends to infinity. But the change is not significant since fi > introduces 
no alteration in our procedure apart from the change of momentum configuration, and it also 
does not interfere with boost-invariance. Therefore, we will limit details of our presentation 
to the simplest massless case, using a momentum configuration that is most convenient when 
fi = 0. We also continue to keep xq in Eq. (3.9b) equal to xiq. For fj, > one has to make 
Xo greater than 0, too. Other choices of the momentum configurations than we adopt here are 
equally possible in the case /i = 0. The one we choose is particularly suitable for executing 
integrals over the loop momenta and extracting the running coupling dependence on the width 
A analytically. Expressions for 'yx{xi, K^g) fo^' ^ are given in Appendix. 

In the massless case, the most convenient configuration is Xio = and K120 = 0^- For massive 
particles one considers 1 > Xio > but there is no compelling reason to consider k^2o 7^ in 0^ 
theory even for massive particles. Note that the parent momentum in Eq. (4.2) is arbitrary 
and not limited by our choice of the momentum configuration (xio,Kf2o)- This feature is not 
readily available when one considers only specific matrix elements of effective Hamiltonians or 
when explicit boost invariance is not preserved by the renormalization group procedure. 

Eq. (4.6) gives 



9x = 9o + 



'24 {Any 



where fx = exp —z^j}^ and /o = exp — z^/Aq. A straightforward integration gives 



(4.8) 



9. = 9.- ^0^(1^ log A , (4.9) 

which exhibits asymptotic freedom. Differentiating with respect to A and keeping terms up to 
order g\ one obtains 

This equation demonstrates the same /? function for coupling constants in effective Hamiltonians 
as obtained in Lagrangian approaches using Feynman diagrams and dimensional regularization, 
when one identifies the renormalization scale with the Hamiltonian width A. This is encour- 
aging but one needs to remember that for comparison of perturbative scattering amplitudes in 
Hamiltonian and Lagrangian approaches it is necessary to make additional calculations and at 
least of fourth order in g. Beyond model matrix studies such as in [|T3], 4th order similarity 



calculations have so far been carried out only in simplified Yukawa model by Maslowski and 
Wi§ckowski ||3^ . 



Integrating Eq. (4.10) one obtains [a = g /A-n 
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^ l + «o(3/327r2)logA/Ao ' ' ' ' 

which shows our result for a boost invariant running couphng constant in effective Hamiltonians. 
Our procedure explains how the running coupling constant can be included in quantum me- 
chanics of effective particles, which is given by the Schrodinger equation with the corresponding 
Hamiltonian Hx. 



5. CONCLUSION 

A remarkable result of the whole procedure in the case of asymptotically free scalar dynamics 
is an extremely simple structure of diverging counterterms and a complete control over involved 
effective interactions, in perturbation theory. Our operator calculus preserves cluster properties 
and allows for evaluation of effective Hamiltonians without limitation to any particular set of 
matrix elements. In other words, we can easily derive integral expressions for matrix elements 
of effective Hamiltonians in the whole Fock space spanned by basis states of effective particles. 
The renormalization group equations are integrated analytically using gaussian similarity form 
factors. Also, the unitary nature of the similarity transformation for effective particle creation 
and annihilation operators removes wave function renormalization constant from the procedure. 
However, the regularization factor introduced in the initial Hamiltonian requires an additional 
finite counterterm that contains a boost invariant function of longitudinal momentum fractions. 

The renormalization group flow for Hamiltonians differs from standard procedures applied 
to S-matrix. The key difference from standard procedures is that the running coupling constant 
is derived as a coefficient in front of a certain term in an effective renormalized light-front Hamil- 
tonian instead of a Lagrangian term or in a scattering amplitude. The advantage of knowing 
an effective Hamiltonian is that not only one can attempt to describe scattering processes but 
also bound states using the corresponding eigenvalue equation. 

Note that the light-front form of relativistic quantum field dynamics re-defines the vacuum 
problem in a way that is only partly understood. [0 But the present calculation shows 
that running coupling constant calculations can be consistently carried out in lowest orders of 
perturbation theory without inclusion of perturbative modifications of the vacuum state. This 
leads to a new demand for similar perturbative analysis of theories that may include effects 
usually thought to be associated with ground state properties, the closest being a scalar theory 
with quartic interaction term and infinitesimally small coupling constant. 

Obviously, a more sophisticated treatment is necessary in gauge theories for many reasons 
but, in particular, because they require an additional cutoff limiting the longitudinal momentum 
fractions from below. The small-x cutoff appears in the unitary transformation U\ as well as 
in the transformations lAxJj'l^, the former depending on and the latter being independent of 
the ultraviolet cutoff A. This difference between the small-x cutoff and the ultraviolet cutoff 
cannot be studied in the case of scalar fields discussed in the present paper. However, the third 
order boost-invariant similarity factors we derived in the case of 0'^ directly apply, for example, 
in calculation of a triple-gluon vertex in effective QCD. 

Finally, we wish to stress the difference between the regularized initial Hamiltonian for 
bare particles and the small width Hamiltonian for effective particles, which contains similarity 
form factors f\. The form factors dampen interactions changing invariant masses by more 
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than A and thus can tame the spread of eigenstate wave functions for low lying eigenvalues 
into regions of high relative momenta of constituents. This feature may lead to exponential 
convergence of the eigenstate expansion in the effective particle basis. Such convergence is 
not expected in the case of bare particles. The fine structure of effective particles would then 
unfold in the transformation UxJA'I^ relating effective degrees of freedom at two different scales, 
one corresponding to the binding scale and the other to the high momentum transfer probe 
in question. This picture encourages opinion that the present calculation provides a pattern 
worth trying in application to realistic theories. 
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Appendix 

The loop integrals in Eq. (2.25) are given in the order implied by first three terms in the 
first bracket. The first term contains the integral /i, the second term in the bracket together 
with the first mass counterterm term contain I2 and the third term in the bracket together 
with the second mass counterterm contain 73. The right-hand side of Eq. (2.25) is given by 
Eq. (3.6) where 7a(^i, k^) — I1 + I2 + I3. In all three terms we have 

Ml2^{l^U+l^')/{XlX2) , (AOOl) 

where 

xi^kf/kt , (A002) 

and 

- xik^ . (A.003) 

The first term integral is 



(AlOl) 



.2/ 2(27r)5 2 [ixi x{l - x){x - Xi) 
where in ^3^(0, c, d) in Eqs. (2.26a-b) one substitutes 

ah = -ha = //^ - Ml^ , (A.102) 

ac = -ca = MI2 - , (A.103) 

ad = -da = MI2 - A*^ , (A.104) 
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with 



and 



be = -cb = MIq - , 
bd — —db — bc/x -\- cd — bajxi + ad , 
cd = -de = Al^ - fjL^ , 

M^^{k^^ + IJ,^)/[x{1-x)] , 

Mie = x''{Ki^^ + li^)/[{x-xi)xi] , 

Mls = xl{4s^ + fi^)/[{x-x,){l-x)] , 

tai = exp |- [?7(a;) + 7]{1 - x)]k-^^ + [r;[(x - Xi)/x] + ?7(a;i/a;)]fi;^6^+ 
+ [77[(X - Xi)/X2] + 77[(1 - x)/x2]]4s^] /A^} , 

i^m = -{I- x)k^^/x2 , 

^bc — ^cb — •^'^3 ' 



p+ _ p+ _ p+ _ p+ _ p+ _ p+ _ 1.+ 



The expression for I2 is 



1 dx 



2] 2(27r)5yo x{l-x) 



-2 



dc X2 



/C2 



d^K' 



rA2 - 2 



1 



4:^3^^(0,6, c, (i) 

+ (1^2) 



^2 

fad-ig ^1^1 



da 2 X2 



{A.105) 
{A.106) 
(A 107) 

(^.108) 
(^.109) 
(yl.llO) 

(Alll) 

(^.112) 
(^.113) 
(^.114) 
(^.115) 



(A.201) 



where in J-'2\o{^j c) in Eq. (2.22) and J^z\{a, b, c, d) in Eqs. (2.26a-b) one substitutes 
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with the same 



and 



ab = -ba ^n'^-M^ , {A.202) 

ac = -co = , (^.203) 

ad = -da = cd , (A.204) 

be = -cb = 6a , (A.205) 

bd = -ci6 = 6c/x2 + cd , (A.206) 

cd = -de = MI2 - , (A.207) 

M^^{k^^ + I^^)/[x{1-x)] , (A208) 

rA2 = cxp {-2[r]{x) + 7]{1 - x)]K^yA^} , (A.209) 

^06 = ^at = = ^2/C3+ , (A210) 

P:d-K, = Pc', = H . (A211) 



The expression for /s is 



/ad ~ 1 

da 



2 

-T+^2Ao(^ C, d) 
/C3 



, /ad 1 r 2 



(A.301) 



where in !F2\o{b, c, d) in Eq. (2.22) and !F^\{a, b, c, d) in Eqs. (2.26a-b) one substitutes 

ab^ad^ MI2 - //^ , (A.302) 

ac = - , (^.303) 

6c = dc = - , (^.304) 

22 



bd = 0, (A.305) 

with the same as for Ji and I2 and 

'"AS = rA2 , (A. 306) 

and all parent momenta equal k^. 
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